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ABSTRACT 

A  method  for  estimating  extinction  coefficients  in  the  near-infrared  and  infrared 
wavelengths  from  satellite  measurements  in  the  visible  was  developed.  Five  tests  were 
devised  to  examine  the  limits  and  sensitivity  of  the  model.  The  first  test  studied  the 
error  inherent  in  the  retrieval  of  the  parameters  which  are  needed  to  describe  the 
distribution  of  atmospheric  particles  and  are  direct  inputs  for  the  calculation  of 
extinction  at  3.75  and  11.00  Jim.  Also  studied  were  errors  associated  with  uncertainties 
in  the  extinction  values,  uncertainties  in  relative  humidity  values,  deviations  of  particle 
size  distribution  from  the  model  and  effects  of  high  winds  on  the  aerosol  distribution. 
Results  indicate  that  the  biggest  error  results  when  wind  generated  aerosols  change  the 
particle  size  distribution  especially  at  radii  larger  than  2  ]im.  The  error  reaches  32% 
for  prediction  at  11.00  fim  at  60%  relative  humidity.  The  smallest  error,  less  than  7% 
for  all  variations,  is  associated  with  the  retrieval  technique  itself.  Errors  up  to  25%  in 
the  measured  satellite  extinction  coefficients  lead  to  errors  of  up  to  25%  in  the 
estimated  values  for  both  the  marine  and  rural  models.  Results  indicate  that  the  rural 
model  at  high  values  of  relative  humidity  is  affected  the  most  with  an  error  of  31%  at  a 
RH  of  95%  at  11.00  |im.  Negative  deviations  in  the  marine  particle  size  distribution 
give  rise  to  large  errors  for  \  =  11.00  um.  For  20%  deviation,  the  error  can  be  as 
high  as  41%.   The  error  decreases  accordingly  as  the  percent  deviation  is  reduced. 
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I.  INTRODUCTION 

A.       MOTIVATION  AND  BACKGROUND 

Electro-optical  surveillance,  guidance  and  weapon  systems  are  affected  by 
oceanographic  and  meteorological  parameters  (Mitchell  et  ai,  1982).  Environmental 
features  such  as  clouds,  ram,  relative  humidity,  haze,  dust,  aerosols  and  precipitation 
have  been  shown  to  have  significant  effects  on  the  performance  of  these  systems. 
Mitchell's  study  points  out  that  systems  such  as  the  REWSON  optical  missile  warning 
set  used  by  helicopters,  medium  range  air-to-surface  missiles,  the  AEGIS  combat 
weapon  system  (a  fully  automatic  surface-to-air  system)  and  high  energy  lasers  have 
been  found  to  be  susceptible  to  the  presence  of  aerosols  in  the  atmosphere.  These 
environmental  parameters  are  not  always  considered  in  evaluating  the  performance  of 
these  weapon  systems.  This  is  because  either  the  environmental  information  is  not 
available,  or  the  scale  of  the  phenomena  is  so  small  that  little  confidence  can  be  placed 
in  the  accuracy  of  its  specification.  In  order  to  improve  the  evaluation  of  the 
performance  of  these  systems  the  applicable  parameters  should  be  included  for  each 
weapon  system. 

Models  such  as  the  LOWTRAN  6  used  for  the  evaluation  of  performance  of 
certain  missile  systems  can  handle  properly  the  absorption  and  scattering  of  radiation 
by  water  vapor  for  the  prediction  of  transmission  losses  but  improvement  could  be 
done  in  the  initialization  to  account  for  a  better  representation  of  aerosols  (Kneizys  et 
ai,  1983).  A  joint  effort  by  the  Naval  Research  Laboratory',  the  Naval  Postgraduate 
School  and  the  Naval  Ocean  Systems  Center  is  presently  underway  to  develop  a  Navy 
Ocean  Vertical  Aerosol  model  (NOVAM)  for  inclusion  into  a  future  version  of 
LOWTRAN.  A  different  approach  involving  small  modifications  of  LOWTRAN  6 
could  be  done  that  would  allow  the  insertion  of  an  extinction  coefficient  value  at  the 
infrared  wavelengths  rather  than  the  actual  parameters  needed  to  calculate  the  particle 
size  distribution.  This  extinction  coefficient  value  would  be  obtained  from  satellite 
estimates  of  optical  depth  and  would  take  into  account  the  aerosols  present  at  the 
target  location  rather  than  at  the  point  from  where  the  weapon  system  is  launched. 
This  is  important  because  some  weapon  systems  rely  on  target  signature  in  determining 
system    performance.      Target     signature     is    related    to     the     target-to-background 


temperature  difference,  which,  when  propagated  through  the  atmosphere,  is  modified 
by  the  transmittance  and  thus  by  the  presence  of  aerosols. 

Since  there  are  weapon  systems  that  operate  at  infrared,  near-infrared  and  visible 
wavelengths,  evaluation  of  the  effect  of  environmental  parameters  must  be  performed 
at  more  than  one  wavelength.  Estimating  the  extinction  of  the  atmosphere  (related  to 
parameters  such  as  visibility)  must  be  done  at  the  visible  and  near-infrared  wavelengths 
and  extrapolated  to  infrared  wavelengths.  This  infrared  extinction  may  be  obtained 
from  satellite  measurements,  but  at  the  infrared  wavelengths  the  power  received  from 
scattered  radiation  is  small  compared  to  visible  wavelengths.  Also,  due  to  the  presence 
of  ambiguities  from  a  variable  background  and  water  vapor,  it  is  more  practical  to 
obtain  measurements  at  visible  wavelengths  and  estimate  the  extinction  at  the  infrared 
wavelengths. 

Studies  by  Gerber  (1985)  show  that  it  is  possible  to  use  visible  light  scattered  by 
aerosols  to  give  useful  estimates  of  extinction  at  infrared  wavelengths.  He  showed  that 
measurements  using  a  combination  of  a  two  channel  integrating  nephelometer.  which 
operates  at  one  wavelength  over  two  angular  ranges,  and  an  axial-scatter  sensor  gave 
the  best  accuracy  ( •**  ±  10%)  for  scaling  over  the  1.06-10.06  Jim  range  for  extinction  by 
marine  aerosols.  However,  his  studies  require  direct  measurement  of  the  scattered  light 
by  the  above  instruments.  The  object  of  this  thesis  is  to  obtain  an  accurate 
representation  of  the  size  distribution  of  the  aerosols  in  the  marine  atmospheric 
boundary  layer  from  satellite  measurements  of  radiance  at  visible  wavelengths.  The 
size  distribution  will  then  be  used  to  obtain  estimates  of  extinction  coefficients  at  near- 
infrared  and  infrared  wavelengths. 

B.       ORGANIZATION 

The  thesis  is  composed  of  five  chapters.  These  include  an  introduction,  theory, 
procedure,  results  and  discussion,  and  a  summary  of  conclusions  and 
recommendations.  Chapter  II  will  discuss  the  theory  of  radiative  scattering  by 
atmospheric  aerosols  and  the  basic  assumptions  used  in  computing  the  scattering  and 
extinction  coefficients.  A  brief  description  of  the  atmospheric  aerosol  models  is  also 
included.  Chapter  III  describes  the  model  for  estimating  extinction  coefficients  at 
near-infrared  and  infrared  wavelengths.  It  discusses  the  assumptions  dealing  with  the 
vertical  distribution  of  atmospheric  particles  in  the  marine  boundary  layer.  Also 
included   is   the   retrieval    of  the    parameters    required   to    describe   the   particle   size 
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distribution  of  the  aerosois  and  the  calculation  of  extinction  at  3.75  and  11.00  Jim. 
Chapter  III  also  discusses  the  validity  of  the  numerical  techniques  used.  Chapter  IV 
discusses  the  procedure  and  results  of  five  tests  developed  to  evaluate  the  sensitivity 
and  limitations  of  the  models.  A  summary  of  the  results  and  recommendations  are 
stated  in  Chapter  V. 
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II.  THEORY 

A.       EXTINCTION  AND  SCATTERING  COEFFICIENTS 

Radiation  can  be  scattered  by  air  molecules,  aerosol  particles  which  are 
suspended  in  the  air  and  by  the  earth's  surface.  Scattering  by  particles  which  are  much 
smaller  than  the  incident  wavelength,  such  as  molecules,  is  described  by  Rayleigh 
theory.  Therefore,  in  the  presence  of  gas  molecules  ( ~  10  ]im)  and  small  particles 
(~0.05  Jim),  radiation  undergoes  Rayleigh  scattering  at  visible  wavelengths.  For 
particles  whose  size  are  comparable  to  or  larger  than  the  wavelength,  the  scattering  is 
referred  to  as  Mie  scattering.  Thus,  in  the  presence  of  cloud  droplets  (~10  Jim) 
radiation  undergoes  Mie  scattering  at  the  infrared  wavelengths,  while  in  the  presence  of 
aerosols  ( —  0. 1- 1.0  Jim)  radiation  undergoes  Mie  scattering  at  the  visible  wavelengths. 
In  a  cloud-free  environment  away  from  large  temperature  or  pressure  gradients  the 
molecular  constituents  do  not  vary  in  the  horizontal  appreciably  and  thus  Rayleigh 
scattering  does  not  vary  significantly  in  the  horizontal.  However,  Mie  scattering 
depends  on  the  size  and  number  of  aerosol  particles.  These  particles  are  primarily 
generated  by  sea  spray  in  a  marine  environment  and  are  composed  primarily  of  dust- 
like particles  over  the  continents.  Advection  of  these  rural  aerosols  would  result  in 
large  horizontal  variations  of  the  particle  size  distribution. 

Scattering  of  light  by  aerosols  is  characterized  by  the  extinction  coefficient.  Pext- 
Pext  is  given  by  the  sum  of  the  extinction  by  scattering  and  absorption. 

Pe.u  =     PScat  +    Pab,  <2-'> 

Absorption  of  radiation  by  aerosols  is  described  by  the  complex  part  of  the  index  of 
refraction.  For  marine  particles  the  complex  part  of  the  index  of  refraction  is  less  than 
10  at  wavelengths  less  than  1  Jim,  and  therefore  absorption  by  marine  aerosols  is  very 
small  (Shettle  and  Fenn,  1979).  For  visible  and  near  infrared  wavelengths,  Pext  for 
marine  aerosols  is  approximately  equal  to  Pscat-  For  higher  wavelengths  and  for  rural 
cases,  the  complex  part  of  the  index  of  refraction  is  of  order  10  so  both  terms  must 
be  included  in  the  calculation  of  extinction  at  3.75  and  11.00  |im.  Gases  in  the 
atmosphere,    such    as    C02    and   water   vapor,    absorb    radiation.     However,    at    the 
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wavelengths  of  interest  (0.63,  0.86.  3.75  and  11.00  Jim)  the  absorption  by  these 
constituents  is  small  and  can  be  neglected  in  our  calculations. 

Given  the  size  and  distribution  of  these  particles,  the  extinction  coefficient  (in 
units  of  per  length)  is  defined  as: 

Pext  =  0l°°7tr2Qext(m'rKdN(r)/dr)dr  »  <2-2> 

where  Qext  is  the  extinction  efficiency  of  a  particle  with  radius  r  and  complex  index  of 
refraction  m,  and  dN/dr  describes  the  size  distribution  of  the  particles.  As  seen  in  Eqn. 
2.2,  the  extinction  coefficient  depends  on  three  terms:  the  cross-sectional  area  (7tr~),  the 
extinction  efficiency  (Qext)  and  the  size  distribution  (dN/dr).  The  relationship  of  r. 
Qext  and  dN/dr  for  a  typical  visual  case  is  presented  in  Fig.  2.1.  Panel  b  of  Fig.  2.1 
shows  a  particle  size  distribution  for  a  relative  humidity  (RH)  value  of  90%  composed 
of  two  modes.  The  modes  correspond  to  small  (mode  1)  and  large  (mode  2)  radii 
particles.  The  nature  of  each  mode  is  described  in  detail  in  IT.  B.  This  example  shows 
that  when  r  is  small,  nr~  and  Qext  are  small,  while  dN/dr  is  large;  when  r  is  large.  kv 
is  large  and  dN/dr  and  Qext  become  small.  The  result  is  that  pext  is  affected  by  a 
bounded  region  of  particle  sizes.  Low  values  of  Qext  bound  p  xt  at  the  small  radii 
while  a  decrease  in  dN/dr  bounds  pext  at  the  large  radii.  For  visible  wavelengths  and  a 
marine  aerosol  size  distribution  model  at  80%  relative  humidity,  ■  studies  done  by 
Durkee  (1987)  include  a  particle  range  from  approximately  0.5  to  5.0  ]im.  To  account 
for  variations  in  wavelength  and  relative  humidity  this  study  incorporates  a  radii  range 
from  0.01  to  50.0  \im  in  the  extinction  coefficient  calculation. 

The  extinction  efficiency  function  is  dependent  on  wavelength  and  relative 
humidity  since  composition  changes  with  increased  condensed  water.  As  wavelength 
increases  Qext  shifts  to  the  right,  weighing  larger  radius  particles  more  than  smaller 
particles  as  depicted  by  Figs.  2.1  and  2.2.  Since  there  are  fewer  particles  at  the  larger 
radii  the  net  effect  of  increasing  wavelength  is  a  decrease  in  the  extinction  coefficient, 
Pexr  Fig.  2.3  shows  the  effect  of  relative  humidity  on  the  size  distribution  of  marine 
particles.  Notice  the  increase  in  the  number  of  large  particles  for  higher  relative 
humidity  values  on  Fig.  2.3.  Since  the  Qexl  function  at  3.75  ]im  weighs  the  larger 
particles  (greater  than  the  mode  radius)  more  than  smaller  particles,  and  there  are  more 
particles  at  the  larger  radii,  the  net  effect  is  an  increase  in  the  cumulative  extinction  for 
higher  relative  humidity  values. 
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Fie.  2.1  a.  Extinction  efficiencv  for  X  =  0.63  Jim  and  RH  =  90%  for  mode  1  (solid) 
antl  mode  2  (dashed),  b.  Model  distribution  of  marine  particles  at  RH  =  90%.  Sum 
mode  1  &  2  (solid),  mode  1  (dashed),  mode  2  (dot),  c.  Cumulative  extinction  from 
Eqn.  2.2. 
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Fie.  2.3    Variation  of  marine  model  size  distribution  of  Shettle  and  Fenn  (1979)  with 
relative  humidity. 
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B.       ATMOSPHERIC  AEROSOL  MODELS 

1.  Introduction 

Atmospheric  aerosol  particles  in  the  atmosphere  vary  greatly  in  their 
concentration,  size  and  composition,  and  consequently  in  their  effects  on  optical  and 
infrared  radiation.  Several  tropospheric  aerosol  models  have  been  developed  by  Shettle 
and  Fenn  (1979)  which  include  the  dependence  of  atmospheric  aerosols  on  relative 
humidity. 

The  size  distribution  for  the  different  aerosol  models  used  here  are  represented 
bv  the  sum  of  two  loe-normal  distributions: 


dN 
dr 


2 
V 


Ni 


i=l\ln(10)   r  <rV27r 


exp 


(log   r  -   log   r . )' 


(2.3) 


where  NT  is  the  cumulative  number  density  of  particles  of  radius  r;  (7  is  the  standard 
deviation;  r.,  NT.  are  the  mode  radius  and  the  total  number  density.  This  form  of  the 
distribution  function  represents  the  multimodal  nature  of  the  atmospheric  aerosols 
which  will  be  discussed  next.  This  study  will  begin  with  a  bimodal  size  distribution 
since  it  has  been  shown  to  be  generally  adequate  to  characterize  the  gross  features  of 
most  aerosol  distributions  (Whitby  and  Cantrell,  1975).  A  third  mode  will  be  added  to 
allow  for  wind  generated,  larger  than  2  urn  mode  radii,  particles.  The  rurai  model  is 
comprised  of  a  small  rural  aerosol  mode  (mode  1)  and  a  large  rural  aerosol  mode 
(mode  2).  The  marine  model  consists  of  a  small  rural  aerosol  mode  (mode  I)  and  an 
oceanic  mode  (mode  2). 

2.   Rural  and  Marine  Aerosol  Models 

The  rural  model  is  intended  to  represent  the  aerosol  under  conditions  where  it 
is  not  directly  influenced  by  urban  and/or  industrial  aerosol  sources.  The  aerosols  are 
assumed  to  be  composed  of  a  mixture  of  70  percent  of  water  soluble  substance  and  30 
percent  dust-like  aerosols.  The  marine  aerosol  model  is  composed  of  a  sea-salt  mode 
and  a  continental  mode  which  was  assumed  to  be  identical  to  the  rural  aerosol  with  the 
exception  that  the  very  large  particles  were  eliminated  since  they  will  be  eventually  lost 
due  to  fallout  as  the  air  masses  move  across  the  oceans.  The  relative  proportions  of. 
aerosols  of  oceanic  or  continental  origins  will  vary,  particularly  in  coastal  regions.    To 


17 


account  for  these  variations,  the  model  permits  adjustment  of  the  relative  contribution 
by  the  oceanic  and  continental  modes.  The  contribution  from  the  marine  size 
distribution  used  in  this  study  varied  from  zero  to  6%  of  the  total  number  of  particles 
of  oceanic  origin  while  the  large  rural  aerosols  varied  from  zero  to  0.025%.  The  small 
radii  particles,  mode  1,  make  up  the  rest  and  were  changed  accordingly. 

As  relative  humidity  increases,  water  vapor  condenses  out  of  the  atmosphere 
onto  the  particulates  suspended  in  the  atmosphere.  This  condensed  water  increases  the 
size  of  the  aerosols  and  changes  their  composition  and  their  effective  refractive  index. 
The  resulting  effect  of  the  aerosols  on  the  absorption  and  scattering  of  light  will  be 
modified  correspondingly.  Table  1  shows  how  the  size  of  a  particle  changes  with 
relative  humidity  for  both  the  rural  and  oceanic  models.  As  relative  humidity 
increases,  the  mode  radii  increases  for  all  modes  and  thus  the  number  of  particles  at  a 
given  radius  increase  with  the  larger  radii  increasing  the  most,  as  seen  previously  in 
Fig.  2.3. 


TABLE 

1 

CHANGES 

IN  AVERAGE 
FUNCTION 

MODE  RADIUS 
OF  RELATIVE 

( rl  AND  r2)  AS  A 
HUMIDITY 

RELATIVE 
HUMIDITY 

RURAL 
rl(jim)    r2(um) 

MARINE 
r2(  urn) 

0% 

0. 02700 

0. 4300 

0. 1600 

50% 

0. 02748 

0. 4377 

0. 1711 

70% 

0. 02846 

0. 4571 

0. 2041 

80% 

0. 03274 

0. 5477 

0. 3180 

90% 

0. 03884 

0. 6462 

0. 3803 

95% 

0. 04238 

0. 7078 

0. 4606 

98% 

0. 04751 

0. 9728 

0. 6024 

99% 

0. 05215 

1.  7555 

0. 7505 

The  marine  model  was  modified  to  allow  the  addition  of  a  wind<  driven  third 
mode  to  the  size  distribution.    As  shown  in  Fig.  2.4  this  third  class  of  aerosol  consists 
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of  the  largest  nuclei  originating  from  the  sea  surface.  It  has  a  mode  radius  of  2  ]xm 
and  the  amplitude  varies  with  wind  speed  (Gathman.  1983).  These  large  particles  can 
be  important  in  the  propagation  of  infrared  radiation  near  the  sea  surface  under  a 
condition  of  high  winds. 


19 


10 
Radius  (um) 


Fig.  2.4  Three-mode  particle  size  distribution  for  marine  model  at  RI1  =  90%  and  A3 
=  1.0.  Sum  of  modes  1  &  2  (solid),  mode  1  (dashed),  mode  2  (dot),  mode  3  (chain- 
dash). 
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III.  PROCEDURE 

A.       LOGIC 

The  main  goal  of  this  study  is  to  be  able  to  estimate  extinction  coefficients  at 
near-infrared  and  infrared  wavelengths  from  satellite  information  at  shorter 
wavelengths.  Channel  1  (0.55  to  .68  Jim)  and  channel  2  (0.70  to  1.00  Jim)  of  the 
NOAA-7  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  return  image 
information  from  which  a  value  of  radiance  is  obtained.  Since  radiance  can  be  shown 
to  be  a  nearly  linear  function  of  aerosol  optical  depth  (Durkee  et  ai,  1986),  the  optical 
depth  due  to  the  presence  of  aerosols  can  be  calculated  as  discussed  in  III.  B. 

The  optical  depth  is  calculated  from  integration  of  the  extinction  coefficient: 

t  =  j'Pextdz  (3.1) 

Above  the  boundary  layer,  as  shown  in  Fig.  3.1  (Fairall  and  Davidson,  1986),  the 
volume  of  aerosols  decreases  due  to  a  decline  in  relative  humidity,  and  therefore  the 
extinction  coefficient,  will  decrease  as  discussed  in  Chapter  II.  A.  Thus,  the  extinction 
above  the  top  of  the  marine  atmospheric  boundary  layer  (MABL)  can  be  considered 
negligible  since  the  number  of  particles  above  this  height  is  significantly  reduced.  The 
integration  is  thus  simplified  and  pext  can  be  approximated  from: 

Pext  =  T/Z  •  <3-2) 

where  Z  is  the  height  of  the  top  of  the  MABL.  Therefore,  Pext  can  be  estimated  from 
satellite  estimates  of  x  and  estimates  of  Z  obtained  from  models  such  as  the  Fleet 
Numerical  Oceanography  Center's  Navy  Operational  Local  Analysis  and  Prediction 
System  (NOLAPS)  (Burk.  and  Thompson,  1982).  Later  in  this  chapter  it  will  be  shown 
that  this  extinction  value  is  the  input  required  in  the  model  in  order  to  retrieve  a 
particle  size  distribution  which  is  used  to  calculate  the  extinction  at  other  wavelengths. 
A  simple  technique  for  estimating  extinction  coefficient  at  longer  wavelengths 
was  attempted  by  plotting  extinction  coefficient  values  in  the  visible  and  near-infrared. 
Fig.  3.2  shows  the  relationship  between  extinction  coefficients  at  0.63  and  3.75  }im  for 
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Fig.  3.1  Height  dependence  of  aerosol  volume  (V)  and  RH  (H)  above  the  ocean  shows 
a  decrease  in  aerosol  above  MABL  due  to  a  decrease  in  RH.  (Fairall  and  Davidson, 
1986). 
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three  values  of  the  contribution  of  the  second  mode,  N2.  Each  distribution  case  is 
varied  by  changing  relative  humidity.  The  total  number  of  particles,  TN,  was  set  at 
4,000  cm"3,  relative  humidity,  RH,  varied  from  0  to  99%,  and  the  contribution  of  the 
second  mode,  N2,  was  set  at  0.0,  0.000125  and  0.000250  for  each  line  respectively.  In 
order  to  find  extinction  values  at  3.75  and  11.00  Jim  an  infinite  number  of  plots  similar 
to  Fig.  3.2  would  have  to  be  generated,  one  for  each  distribution  shape.  A  simple 
prediction  using  one  wavelength  is  therefore  impossible.  This  result  is  consistent  with 
Gerber  (1985)  who  showed  for  various  atmospheric  simulations  that  no  correlation 
exists  between  extinction  at  visible  and  near-infra  red  wavelengths  when  scaling  with 
only  one  wavelength. 

We  have  shown  that  the  extinction  value  depends  on  the  number  of  particles 
present  and  on  their  distribution.  Thus,  a  method  using  satellite-detected  radiance  at 
two  wavelengths  is  proposed  for  estimating  extinction  at  the  3.75  and  11.00  Jim 
wavelengths.  Approximations  made  in  order  to  convert  radiance  measurements  from 
satellite  data  to  an  extinction  value  are  discussed  in  III.  B.  From  an  input  of 
extinction  values  at  0.63  and  0.86  ^m  the  model  can  retrieve  the  total  number  of 
particles  (TN)  and  the  percentage  of  the  second  mode  (N2).  With  knowledge  of  these 
parameters,  which  are  needed  to  describe  the  particle  size  distribution,  and  an  effective 
humidity  value  for  the  boundary  layer,  the  model  can  calculate  extinction  at  3.75  and 
11.00  fim  as  discussed  in  III.  C. 

B.       CALCULATION  OF  OPTICAL  DEPTH  AND  EXTINCTION  COEFFICIENT 

Channel  1  (0.63  Jim)  and  channel  2  (0.86  Jim)  of  the  NOAA-7  AVHRR  satellite 
returns  image  information  in  the  form  of  digital  brightness  counts.  A  digital  number 
(which  is  proportional  to  the  amount  of  radiation  received  at  a  satellite)  is  assigned  to 
each  image  element.  This  inserts  some  uncertainty  since  only  integers  are  used  in  the 
digitizing  process.  Fig.  3.3  shows  the  conversion  between  digital  counts  and  physical 
units  for  the  AVHRR  (Durkee,  1984). 

The  radiance  at  satellite  altitude  (L«j)  consists  of  contributions  from  the  following 
terms: 

Ls=  (Lw  +  LG)T  +  LR  +  LA,  (3.3) 
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Fie.   3.2    Estimation  of  extinction  coefficient  at   3.75  Jim  from  0.63   [im  for  three 
different  particle  size  distribution,  described  by  N2. 
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Fie.  3.3    Conversion  of  digital  counts  to  radiance  for  the  AVHRR  sensors  (Durkee, 
19S4).  & 
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where  L^y  is  the  water  leaving'  radiance  caused  by  subsurface  reflectance,  Lq  is  the 
contribution  from  the  specular  reflection  off  the  sea  surface,  L^  is  the  path-added 
radiance  due  to  Rayleigh  scattering  in  the  atmosphere,  L^  is  the  path-added  radiance 
due  to  Mie  scattering  in  the  marine  boundary  layer  and  T  is  the  atmospheric 
transmittance. 

For  wavelengths  in  the  red-visible  to  near-infrared,  L^y  generally  will  be  a  small 
contribution  to  upwelling  radiance.  Lq,  also  referred  to  as  the  radiance  added  by 
'sunglint',  can  be  estimated  from  knowledge  of  the  sun-earth-satellite  geometry  and 
surface  roughness.  Lq  is  negligible  under  all  but  certain  geometries.  L^  is  obtained 
from  knowledge  of  the  Rayleigh  optical  depth  (which  can  be  estimated  from  profiles  of 
pressure  and  temperature)  and  scene  geometry  (McCartney,  1976). 

In  the  single  scattering  approximation  we  consider  the  radiation  which  is 
scattered  only  once  by  marine  aerosols.  If  we  assume  that  the  upward  intensity  at  the 
bottom  of  the  atmosphere  is  zero  and  for  atmospheres  with  small  optical  depth,  it  can 
be  shown  that  radiance  is  a  linear  function  of  optical  depth  (Durkee  el  a/.,  1986).  The 
optical  depth  due  to  the  presence  of  aerosols  can  then  be  estimated  from  the  aerosol 
radiance.  Eqn.  3.2  can  then  be  used  to  calculate  the  extinction  coefficient  required  for 
input  into  the  model.  For  the  purpose  of  this  study  it  is  assumed  that  the  extinction 
values  at  0.63  and  0.86  fim  are  available  from  estimates  of  optical  depth  (T)  and 
boundary  layer  depth  (Z). 

C.       ESTIMATION  OF  IR  EXTINCTION  COEFFICIENT 

Fig.  3.4  shows  a  basic  outline  of  the  model  for  estimating  extinction  at  near- 
infrared  and  infrared  wavelengths.  This  model  can  be  divided  into  two  sections.  In  the 
first  part  of  the  model,  input  of  the  aerosol  model  type  (marine  or  rural)  and  the 
effective  relative  humidity  for  the  MABL  defines  the  location  of  the  average  mode 
radius,  rl  and  r2.  Input  of  the  extinction  coefficient  at  0.63  and  0.S6  Jim  will  determine 
the  magnitude  of  the  retrieved  total  number  of  particles,  TN,  and  percentage 
distribution  of  the  second  mode,  N2.  These  are  the  parameters  necessary  to  describe 
the  distribution  of  the  aerosols  particles.  The  second  part  of  the  model  uses  these 
values  of  TN  and  N2  and  an  effective  relative  humidity  value  for  the  boundary  layer  to 
estimate  the  extinction  coefficient  at  3.75  and  11.00  [lm. 
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Fig.    3.4     Outline    for    the    estimation    of  extinction    coefficients    at    NIR    and    IR 
wavelengths. 
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1 .  Retrieval  of  TN  ana  N2 

This  portion  of  the  mode!  is  used  to  calculate  the  total  number  of  particles 
and  the  contribution  of  the  second  mode  needed  to  describe  the  aerosol  distribution 
using  a  bilinear  interpolation.  The  model's  application  is  limited  by  the  total  number 
of  particles  and  the  variations  allowed.  For  this  study,  in  the  marine  case,  TN  was 
allowed  to  vary  from  1,000  to  7,000  cm"J  and  the  second  mode  contribution,  N2, 
varied  from  0.0  co  0.06.  For  the  rural  case.  TN  was  allowed  to  vary  from  10,000  to 
20,000  cm  and  the  second  mode  contribution,  N2,  varied  from  0.0  to  0.000250. 
These  are  the  current  limitations  of  the  model  but  the  model  can  be  easily  expanded  to 
allow  for  greater  variations.  The  data  files  are  comprised  of  7  X  7  Pext  matrices  for 
marine  and  rural  cases  at  0.63  and  0.86  ]im,  and  for  RH  =  0,  50,  70,  80,  90,  95,  98  and 
99%.  These  matrices  were  obtained  by  running  the  second  portion  of  the  model  for 
the  following  values  of  TN  and  N2.  For  the  marine  case,  N2  values  used  were  0.0, 
0.01,  0.02,  0.03,  0.04,  0.05  and  0.06.  For  the  rural  case  N2  values  used  were  0.0. 
0.000050,  0.000100,  0.000125.  0.000150.  0.000200,  0.000250.  In  order  to  run  the  model. 
an  effective  humidity  for  the  boundary  layer  and  values  of  extinction  at  0.63  and  0.86 
\im  must  be  entered.  A  subroutine  in  the  model,  using  a  Hermite  interpolation,  allows 
calculation  of  a  7  X  7  pext  matrix  for  any  integer  value  of  relative  humidity  less  than 
100%  not  in  the  data  files.  .    . 

As  discussed  in  II.  A..  Pext  is  a  function  of  TN,  N2,  RH  and  \.  An  infinite 
number  of  contributions  of  TN  and  N2  will  give  rise  to  a  predetermined  extinction 
coefficient  at  a  given  wavelength.  But  only  one  TN  and  N2  combination  will  give  rise 
to  a  given  Pext  at  X1  and  a  given  Pext  at  "k2  Therefore,  we  must  generate  lines 
describing  these  combinations  at  each  wavelength  and  solve  for  the  intersection  of 
these  two  lines  which  will  give  the  value  of  TN  and  N2  associated  with  the  input  pext 
at  the  two  visible  wavelengths. 

The  following  description  shows  how  the  model  mathematically  solves  for  the 
two  parameters,  TN  and  N2,  needed  to  describe  the  aerosol  distribution  in  question. 
For  a  given  extinction  value,  a  corresponding  TN  value  is  calculated  at  each  N2.  This 
generates  a  set  of  seven  values  of  TN,  one  at  each  N2  value,  for  each  X.  The 
intersection  between  the  curves  formed  by  these  sets  (Figs.  3.5  and  3.6)  represents  the 
best  approximation  of  TN  and  N2  which  describes  the  particle  distribution 
corresponding  to  the  extinction  coefficient  value  observed  at  0.63  and  0.86  |im..  This 
single   value   of  TN  and  N2  is   obtained   by  a  bilinear  interpolation.     The   bilinear 
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interpolation  consists  of  obtaining  an  equation  for  each  line  using  a  point-form  for  the 
nearest  point  on  each  side  of  the  intersection.  By  solving  the  two  equations  the 
intersection  of  the  two  curves  is  determined  and  therefore  an  estimate  of  TN  and  N2  is 
obtained. 

2.  Calculation  of  Pext  at  3.75  and  11.00  Jim 

This  portion  of  the  model  calculates  an  extinction  coefficient  value,  Pext,  at 
near-infrared  and  infrared  wavelengths,  from  knowledge  of  the  particle  size  distribution 
and  an  effective  value  of  relative  humidity  for  the  marine  boundary  layer.  Data  files 
containing  mode  radii  for  marine  and  rural  cases  are  supplied.  These  values  have  been 
interpolated  using  a  Hermite  interpolation  for  integer  values  of  relative  humidity 
ranging  from  0  to  99%.  Also  included  are  data  files  containing  values  of  index  of 
refraction  for  marine  and  rural  cases.  These  values  have  also  been  interpolated  using  a 
Hermite  interpolation  for  integer  values  of  relative  humidity  ranging  from  0  to  99%. 
Uninterpolated  values  for  mode  radii  and  index  of  refraction  at  RH  =  0,  50,  70,  SO.  90, 
95,  96.  97,  98  and  99%  were  obtained  from  Shettle  and  Fenn  (1979).  The  value  of  TN 
and  N2  obtained  from  the  previous  portion  of  the  model  and  an  effective  relative 
humidity,  RH,  for  the  marine  boundary  layer  are  entered  into  the  model.  Using  Eqns. 
2.3  and  2.2  the  model  calculates  (3     t  at  3.75  and  1 1.00  urn. 

D.       VALIDITY  OF  NUMERICAL  TECHNIQUES 

Values  for  the  average  radius  of  each  mode,  rl  and  r2.  and  the  real  and 
imaginary  parts  of  the  index  of  refraction  for  the  marine  and  rural  cases  are  provided 
for  only  certain  relative  humidity  values  as  discussed  in  the  previous  section.  Figs.  3.7 
and  3.S  are  a  plot  of  the  average  mode  radius  as  a  function  of  relative  humidity  and 
show  a  monotonic  behavior  for  both  rl  and  r2.  Due  to  the  flat  section  followed  by  a 
rapid  nonlinear  portion  for  relative  humidity  values  greater  than  70%,  a  Hermite  cubic 
numerical  interpolation  was  used  to  generate  values  of  rl  and  r2  for  relative  humidity 
values  ranging  from  0  to  99%.  Figures  3.9  and  3.10  show  a  similar  monotonic  trend 
for  the  scattering  extinction  for  rural  and  marine  cases  as  a  function  of  relative 
humidity.  Based  on  this  result,  matrices  needed  for  extinction  values  during  the 
retrieval  portion  of  the  model  are  generated  by  a  Hermitic  interpolation.  The  real  and 
imaginary  portions  of  the  index  of  refraction  also  follow  this  behavior  and  a  Hermitic 
interpolation  was  therefore  employed.  The  Hermitic  cubic  function  in  this  program 
was  developed  by  the  Lawrence  Livermore  National  Laboratory  (Fritsch,  1982).    Such 


29 


7.0 

e.a- 

6.2- 

5.8- 

O 

o 

d      5.4- 

X 


c 

c 
o 


5.0 


4.5- 


<a      4.2- 


o 

O        3.8- 

C 

o 

U        3.4- 
<D 

~    3.0- 

a 

a. 

"3     2fl" 


2.2 

1.8- 
1.4- 


1.0 


TN  V8  N2  AT  0.83  Mil 
TN  V8  N2  AT  0.88  urn 


1 1 1 1 1 1 1 1  1 1 1 

0.0     0.5       1.0       1.5       2.0       2.5       3.0      3.5      4.0      4.5      6.0      6.6      6.0 

Percent  Mode  2  x  100 


Fie.  3.5  TX  value  for  corresponding  N2  at  0.63  {solid)  and  0.S6  (dashed)  Jim 
wavelength.  Intersection  gives  best  estimate  of  IN  and  N2  lor  aerosol  marine 
distribution  given  by  input  extinction  coellicients. 


30 


2.00 

1.95 

1.90 

1.86 

r- 

1.80 

o 

o 

o 

1.75 

o 

X 

1.70 

— » 

E 

1.65 

cj 

%~' 

1.60 

c 

o 

*- 

1.65 

CO 

k_ 

-*- 

c 

1.60 

o 

o 

c 

1.45 

o 

o 

1.40 

© 

o 

1.36 

**> 

km 

1.30 

CO 

1.26 

.0 

h- 

1.20 

1.16 

1.10 

1.05 

1.00 

TN  V8  N2  AT  0.83  ^m 
TN  V8  N2  AT  0.88  pm 


0.00    0.26      0.50 


— I 1 

0.75       1.00 


1.25       1.50       1.75       2.00      2.25       2.60 


Percent  Mode  2  x  10000 


Fig.  3.6  TN  value  for  corresponding  N2  at  0.63  (solid)  and  0.86  (dashed)  Jim 
wavelength.  Intersection  cives  best  estimate  of  IN  and  N2  for  the  rural  aerosol 
distribution  given  by  input  extinction  cocdicients. 


31 


an  interpolating  scheme  is  more  reasonable  than  a  cubic  spline  if  the  data  contain  both 
"steep"  and  "flat"  sections  as  shown  in  Figs.  3.7  -  3.10. 

As  discussed  in  the  previous  section,  a  bilinear  interpolation  was  used  to  solve  for 
the  intersection  of  the  two  lines  and  retrieve  the  parameters,  TN  and  N2,  needed  to 
describe  the  particle  size  distribution  of  aerosols.  This  technique  was  appropriate  for 
this  case,  as  seen  in  Figs.  3.5  and  3.6,  since  the  equation  of  the  lines  are  generated  with 
points  in  the  immediate  vicinity  of  the  intersection  and  a  linear  approximation  between 
such  short  distances  can  be  made.  Fig.  3.5  shows  that  if  the  points  were  separated  a 
greater  distance  a  linear  approximation  would  not  hold  true  for  the  marine  case  since 
over  the  interval  used  for  N2  (0.0  to  0.06)  the  curve  is  not  linear.  However,  it  would 
still  hold  true  for  the  rural  case  since  the  lines  are  almost  linear  as  seen  in  Fig.  3.6.  For 
the  rural  case  the  values  immediately  following  the  intersection  (A1-A4)  were  used  to 
find  the  intersection  point.  A  set  of  values  (B1-B4)  further  away  from  the  intersection 
were  used  and  compared  with  the  first  set.  Prediction  errors  for  extinction  at  3.75  and 
11.00  Jim  were  calculated.  The  results  shown  in  Table  2  show  excellent  agreement 
since  the  errors  are  less  than  1%  for  both  infrared  wavelengths.  Thus  the  bilinear 
interpolation  is  suitable  for  the  rural  case  and  in  this  study  an  acceptable  method  of 
interpolation  for  the  marine  case. 


TABLE  2 

BILINEAR  INTERPOLATION  VALIDITY  TEST 

pextfkm"1) 

ERROR(%) 

3. 75^m 

11. OOum 

3. 75^m      11. OOum 

INPUT 

TN    15575cm"3 
N2    0. 000115 

0. 02446 

0. 01595 

-- 

A 

TN   15551cm"3 
N2    0. 000114 

0. 02443 

0. 01593 

-0.  12         -0. 12 

B 

TN    15550cm'3 
N2    0.000112 

0. 02443 

0.01593 

-0.  12         -0.  12 
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models  from  Table  1. 
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IV.  RESULTS  AND  DISCUSSION 

After  developing  a  technique  to  estimate  extinction  values  in  the  infrared 
wavelengths  from  values  in  the  visible  and  near-infrared,  it  is  necessary  to  t.est  its  limits 
and  sensitivity  to  several  parameters  that  are  included  in  the  model.  Five  tests  were 
designed  with  this  purpose.  The  first  test  is  a  basic  analysis  of  the  technique  itself  in 
the  retrieval  of  the  particle  size  distribution.  Two  of  the  tests  study  how  variations  in 
the  input  parameters,  extinction  coefficient  and  relative  humidity  affect  the  estimated 
value.  The  fourth  test  looks  at  how  deviations  from  the  model  size  distributions  affect 
the  estimated  extinction  values.  The  last  test  compares  how  an  extinction  obtained 
from  a  three-mode  particle  size  distribution  compares  to  an  extinction  calculated  from 
a  model  which  assumes  a  two-mode  particle  size  distribution. 

A.       BASIC  RETRIEVAL  ANALYSIS 

The  first  test  examines  the  error  associated  with  the  retrieval  of  a  "perfect" 
satellite  estimate  of  extinction  after  use  of  the  interpolating  technique.  The  procedure 
consists  of  the  following  steps: 

1.  Calculate  an  initial  extinction  values  at  0.63  and  0.86  urn.  for  a  given  TN,  N2 
and  RH. 

2.  Input  extinction  values  into  the  first  portion  of  the  model  to  retrieve  TN  and 
N2  which  describe  the  particle  size  distribution. 

3.  Input  TN  and  N2  into  the  second  section  of  the  model  to  obtain  an  estimate  of 
the  extinction  values  for  all  four  wavelengths  of  interest:  0.63,  0.86.  3.75  and 
11.00  ^im 

4.  Calculate  the  error  from: 

error  =   100(J3ext  estimated  -  Pext  initial);  (3ext  initial  (4.1) 

5.  The  procedure  is  done  for  both  the  marine  and  rural  cases. 

The  basic  retrieval  test  for  the  marine  case  was  run  with  values  of  TN  and  N2  set 
to  3,500  cm  and  0.035  respectively.  In  the  rural  case  TN  was  set  to  12,000  cm°  and 
N2  to  0.000175.  Humidity  values  of  25,  50,  75,  85,  96,  98  and  99%  were  used. 
Relative  humidity  values  of  25,  75  and  85%  needed  interpolation  by  the  model.  As 
seen  from  the  results  in  Table  3  the  error  associated  with  the  retrieval  technique  is  less 
than  2.5%  for  relative  humidity  values  less  than  99%.    However,  this  value  increases  as 
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me  number  of  parameters  that  need  to  be  interpolated  increases.  With  all  parameters 
interpolated,  the  error  is  still  less  than  10%.  A  larger  error  at  the  higher  humidities 
was  expected  because  this  is  where  extinction  changes  the  most  as  shown  in  Figs.  3.9 
and  3.10.  The  larger  error  at  a  relative  humidity  value  of  25%  is  probably  due  to  the 
fac:  that  the  interpolation  for  this  value  occurs  over  a  large  range,  between  0  and  50% 
whereas  at  the  higher  humidities  the  difference  is  10%  or  less.  The  error  associated 
with  a  relative  humidity  value  of  50%  was  the  lowest  and  uses  an  uninterpolated 
matrix  at  a  relative  humidity  where  there  is  very  little  variation  in  extinction.  This 
shows  that  the  error  due  to  interpolation  could  be  lowered  by  decreasing  the  range 
over  which  the  interpolation  is  to  be  performed.  This  can  be  accomplished  by 
expanding  the  number  of  relative  humidity  values  that  are  used  by  the  model. 


TABLE  3 

ERROR  ASSOCIATED  WITH  BASIC  RETRIEVAL  TECHNIQUE 

a.    MARINE 

MODEL 

RH(%) 

RETRIEVAL 

ERROR(%) 

3. 75^m 

11.  OO^im 

25 
75 
96 
98 

1.  9 
-0.  2 

0.  7 

1.  9 

2.  1 
-0.  2 

0.  7 

1.  9 

b.    RURAL 

MODEL 

RH(%) 

RETRIEVAL 

ERROR(%) 

3. 75nm 

11.  Q0]im 

0 
50 
85 
99 

0.  5 

0.  1 

1.  8 
6.5 

0.  3 
0.  7 
1.5 
6.  5 

B.       EXTINCTION  ERROR  ANALYSIS 

The  second  test  evaluated  the  sensitivity  of  the  technique  to  errors  in  the 
extinction  values  used  as  inputs.  Several  assumptions  were  made,  as  discussed  in 
Chapter  II,   for  obtaining   the  initial  input  extinction   value  from  satellite  radiance 
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measurements.  Therefore,  the  extinction  coefficient  estimates  may  have  significant 
errors.  The  test  consisted  of  assuming  the  initial  extinction  value  is  in  error  by  ±5,  10. 
20%.   The  procedure  was  as  follows: 

1.  For  a  given  TN,  N2,  and  RH,  calculate  Pext  at  the  four  wavelengths  of  interest. 

2.  Using  the  value  obtained  in  step  (1)  calculate  a  B.Yt  value  at  0.63  and  0.86  urn 
that  is  off  by  ±5%.  ext 

3.  Input  above  value  into  the  first  section  of  the  model  to  obtain  a  retrieved  TN 
and  N2. 

4.  Input  the  retrieved  TN  and  N2  value  into  the  second  portion  of  the  model  to 
obtain  an  estimated  extinction  value  at  3.75  and  11.00  ^im. 

5.  Calculate  the  error  as  discussed  in  IV.  A. 

6.  Repeat  procedure  for  ±10  and  ±20%. 

The  marine  model  was  run  for  a  TN  value  of  3,500  cm"3  and  an  N2  value  of 
0.035.  For  the  rural  case  TN  was  set  at  15,500  cm  and  N2  at  0.000135,  using  relative 
humidities  of  25,  75,  96  and  98%.  The  extinction  values  were  then  varied  by  ±5.  10, 
20  and  25%,  as  explained  above  for  each  humidity  value.  Table  4  shows  the  error 
associated  with  a  variation  in  RH  on  an  estimated  extinction  coefficient  value  at  near- 
infrared  and  infrared  wavelengths.  As  seen  from  the  results,  the  magnitude  of  the  error 
is  basically  determined  from  the  accuracy  of  the  extinction  value  within  a  couple  of 
percentage  points.  This  was  expected  since  the  retrieval  method  itself  is  accurate  to 
better  than  2%  for  most  relative  humidities.  Slightly  larger  errors  (3-5%)  for  the  rural 
model  are  observed  for  a  humidity  value  of  98%  where  the  value  of  extinction  grows 
exponentially  as  seen  in  Fig.  3.10.  The  larger  errors  for  the  rural  model  at  RH  =  98% 
are  observed  for  extinction  errors  of  20%  or  larger.  This  follows  the  inherent  error  of 
the  retrieval  technique  at  this  humidity.  Thus  error  in  the  predicted  extinction  is  less 
than  3%  for  all  humidities  for  both  models  after  considering  the  initial  extinction  error 
and  the  error  inherent  in  the  retrieval  technique. 

C.       RELATIVE  HUMIDITY  ERROR  ANALYSIS 

Several  assumptions  were  made  in  the  development  of  the  model  for  estimating 
extinction  coefficients  relating  to  relative  humidity.  The  first  assumption  deals  with  the 
fact  that  relative  humidity  varies  in  the  marine  atmospheric  boundary  layer  (MABL) 
with  height.  Relative  humidity  generally  increases  from  the  surface  to  the  top  of  the 
MABL  then  decreases  drastically  to  near  zero  at  the  top  of  the  inversion  layer  as 
previously  shown  in  Fig.  3.1.  This  model  assumes  input  of  an  average  value  of  relative 
humidity.   The  second  assumption  is  that  the  relative  humidity  value  is  known.    It  may 
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TABLE  4 

ERROR  ASSOCIATED  WITH  UNCERTAINTY  IN  EXTINCTION 

VALUE 


MARINE   MODEL 


RH(%)    %ERROR   ON   pext 


PREDICTED   ERROR(%) 
3.  75\im  11.  OO^im 


25 
25 
25 
25 

±    5 
±10 
±20 
±25 

7.  1 
11.  8 
22.  3 

27.  4 

-    3.  2 
-8.3 
-18.  4 
-23.  7 

6.  7 
11.  9 
22.  2 
27.  3 

-3.1 
-  8.  3 
-18.  6 
-23.  7 

75 
75 
75 
75 

±    5 
±10 
±20 
±25 

4.  8 

9.  8 

19.  5 

24.  5 

-    5.  3 
-10.  2 
-20.  1 
-25.  1 

4.  9 

9.  8 

19.  8 

24.  7 

-  5.  3 
-10.  2 
-20.  2 
-25.  0 

96 
96 
96 
96 

±    5 
±10 
±20 
±25 

5.  7 
10.  7 

-    4.  4 
-9.4 
-19.  5 
-24.  5 

5.  8 
10.  8 

-  4.  3 

-  9.  3 
-19.  4 
-24.  3 

98 
98 
98 

98 

±    5 
±10 
±20 

±25 

7.  0 
12.  1 
22.  3 
27.  4 

-  3.  2 

-  8.  3 
-18.  5 
-23.  7 

7.  0 
12.  1 
22.  3 
27.  4 

-  3.  2 

-  8.  3 
-18.  5 
-23.  6 

b.  RURAL  MODEL 


RH(%)  %ERR0R  ON  pext 


PREDICTED  ERR0R(%) 


25 

25 
25 
25 

75 
75 
75 
75 

96 
96 
96 
96 

98 
98 
98 
98 


±  5 

±10 
±20 
±25 

±  5 
±10 
±20 
±25 

±  5 
±10 
±20 
±25 

±  5 
±10 
±20 
±25 


3.  75]im 


5.  4 

-    5. 

-i 

11.  5 

-10. 

9 

19.  4 

-18. 

8 

24.  3 

-24. 

2 

5.  3 

-    4. 

9 

10.2 

-    9. 

7 

20.  4 

-19. 

9 

25.2 

-25. 

2 

3.  3 

-    8. 

9 

8.  4 

-13. 

6 

18.2 

-23. 

6 

23.  1 

-28. 

3 

9.  7 

-    7. 

8 

10.  2 

-    6. 

5 

25.  4 

-16. 

8 

29.  9 

-22. 

4 

11. OOum 


5.  8 

11.  5 

19.  4 
25.  2 

5.  6 
10.  5 

20.  4 
25.9 

3.  2 

8.  1 
18.  0 
22.  9 

9.  5 
14.  9 
25.  3 
29.  7 


-  5.  8 

-10.  8 
-19.  4 
-24.  5 

-  4.  3 

-  9.  3 
-19.  8 
-24.  7 

-  8.  7 
-13.  6 
-23.  5 
-28.  4 

-  7.  3 

-  6.  6 
-16.  9 
-22.  4 
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be__  obtained  from  in  situ  measurements  or  in  the  future,  perhaps,  from  satellite  data 
(Kren,  1987).  Changes  in  the  relative  humidity  change  the  shape  of  the  particle  size 
distribution  as  shown  in  Fig.  2.3.  In  this  study  it  was  assumed  that  estimates  of 
relative  humidity  value  can  be  obtained  within  a  certain  degree  of  accuracy  ( ±5%).  A 
series  of  commercial  humidity  sensors  were  tested  under  laboratory  conditions  (Muller 
and  Beekman,  1987).  Test  results  show  that  for  some  of  the  sensors  accuracies  of  95% 
or  better  were  achieved.  Therefore,  the  test  performed  was  in  keeping  with  the 
accuracy  reported  for  humidity  sensors. 

The  procedure  for  this  study  consisted  of  the  following: 

1.  For  a  given  TN,  N2  and  RH.  calculate  the  initial  extinction;  for  example  enter 
an  'observed'  (the  real  RH  value)  RH  =  65%. 

2.  Use  the  extinction  values  obtained  for  0.63  and  0.S6  \im  and  input  into  the  first 
section  of  the  model.  Enter  an  assumed'  RH  value  that  is  +o%  (for  example 
enter  RH  =  70%)  to  obtain  TN  and  \"2. 

3.  Enter  TN  and  N2  from  above  and  the  same  RH  value  used  in  step  2  into  the 
second  portion  of  the  model  to  obtain  the  estimated  extinction.  Followma  our 
example,  enter  an  'assumed'  (due  to  errors  in  the  instrument  that  measures~RH) 
RH  =  70%. 

4.  Calculate  error  as  discussed  in  IV.  A. 

5.  Repeat  for  -5%.  Also  repeat  for  different  RH  values  for  both  marine  and  rural 
models. 

The  marine  case  was  run  for  TN  set  to  3.500  cm'J  and  N2  set  to  0.0175  whereas 

the  rural  case  values  were  15.000. cm" J  and  0.000135  respectively.    Table  5  shows  that. 

for  the  marine  model,  a  5%  error  in  relative  humidity  values  give  rise  to  a  range  of 

5-14%  errors  in  J3ext  at  the  near-infrared  and  infrared  wavelengths.    In  the  rural  model 

a  similar  error  in  RH  gives  rise  to  estimated  J3ext  values  being  in  error  from  9-31%.   As 

seen  in  Table  5,  the  rural  model  is  more  sensitive  to  errors  in  relative  humidity  at  11.00 

pm  wavelength   than  the   marine   model.     Recall  that   extinction  coefficients   are  a 

function  of  the  particle  size  distribution.   At  0.63  and  0.86  Jim  the  satellite  data  contain 

more  information  regarding  the  first  mode  in  the  rural  case  in  contrast  to  more 

information  on  the  second  mode  for  the  marine  case.   This  can  be  seen  if  a  vertical  line 

is  drawn  for  radii  of  0.63  and  0.86  fim  on  Figs.  4.1  and  4.2  until  they  intersect  the  solid 

line.    The  slope  of  the  line  at  the  point  of  intersection  for  the  marine  case  is  more 

dependent  on  changes  in  the  shape  of  the  second  mode.    Thus,  a  decrease  in  the  slope 

infers  a  larger  number  of  particles  at  the  larger  radii.   As  relative  humidity  increases  the 

number  of  particles  at  the  larger  radii  increases.   At  the  higher  wavelengths  these  larger 

particles  contribute  more  to   the  extinction  coefficient   since  they  are  more  heavily 
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weighted.    Since  the  information  in  Pext  about  the  second  mode  is  limited  in  the  rural 
model,  a  larger  error  is  expected  as  observed  in  Table  5. 

D.       PARTICLE  SIZE  DISTRIBUTION  ERROR  ANALYSIS 

The  next  study  investigates  the  error  associated  with  errors  in  the  particle  size 
distribution  obtained  by  changing  the  average  radius  of  the  second  mode.  The  particle 
size  distribution  can  be  altered  by  changes  in  the  relative  humidity,  as  was  shown  in 
III.  C.  The  distribution  can  also  be  changed  by  changing  the  average  radius  of  the 
particles  for  the  second  mode.  r2,  of  Eqn.  2.3.  The  second  mode  of  the  marine  case 
was  the  only  one  tested  since  it  is  known  that  it  consists  of  sea  spray  particles,  which 
are  more  likely  to  vary7  in  size  due  to  variations  in  the  wind.  Changing  r2  by  ±20%, 
as  seen  on  Fig.  4.3,  moves  the  center  of  the  second  mode.  The  model  was  run  with  TN 
set  to  4.000  cm"J  and  N2  set  to  0.02.   The  procedure  is  as  follows: 

1.  Chanee  r2  value  for  a  aiven  relative  humidity  value  in  the  prediction  model  by 
±5%.  Calculate  the  initial  extinction  values  representing  a  non-ideal  particle 
size  distribution.    Reset  the  parameters. 

2.  Enter  extinction  values  into  the  first  portion  of  the  model  to  obtain  values  for 
TN  and  N2. 

3.  Enter  TN  and  N2  into  the  second  section  of  the  model  to  obtain  the  predicted 
extinction  values. 

4.  Calculate  error  as  discussed  in  IV.  A. 

5.  Repeat  for  ±10%  and  ±20%  r2  variations  in  the  marine  case  only. 

As  we  increase  r2  we-  change  the  distribution  by  increasing  the  number  of 
particles  with  a  larger  radii.  Table  6  shows  that  as  we  increase  the  variation  in  r2  the 
error  associated  with  the  positive  deviations  at  X  =  3.75  Jim  increases  from  0.5  to  10% 
while  the  error  associated  with  the  negative  deviations  increased  from  6  to  25%.  For 
X=  11.00  urn  the  error  associated  with  the  positive  deviations  increased  from  0  to  20% 
whereas  for  negative  deviations  the  error  increased  from  9  to  41%.  The  error  is 
greatest  for  negative  deviations  at  the  11.00  ^m  wavelength.  Again,  at  the  0.63  and 
0.86  ]im  wavelengths  the  extinction  efficiency  function  weighs  the  larger  particles  more 
heavily  for  higher  wavelengths  as  shown  in  Figs.  2.1  and  2.2.  If  a  vertical  line  is  drawn 
to  the  solid  line  for  X  =  0.63  and  0.86  um,  as  seen  in  Fig.  4.3,  it  can  be  shown  that 
more  information  is  contained  in  Pext  at  0.63  and  0.86  Jim  about  the  second  mode  for 
positive  deviations.   The  net  result  is  that  the  error  is  larger  for  the  -20%  r2  change. 
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Fig.  4.1  Effect  on  the  distribution  of  marine  particles  by  changing  RH  by  20%. 
Vertical  lines  at  Qext  peak  for  A.  =  0.63  and  0.S.6  Jim.  Sum  of  modes  1  &  2  (solid), 
mode  1  (dashed),  mode  2  (dot). 
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Fig.  4.2  Effect  on  the  distribution  of  rural  particles  by  changing  RH  by  20%.  Vertical 
lines  at  Qext  peak  for  \  =  0.63  and  0.S6  fim.  Sum  of  modes  l  &  2  (solid),  mode  1 
(dashed),  mode  2  (dot). 
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TABLE  5 

RELATIVE  HUMIDITY  SENSITIVITY  STUDY 

a.    MARINE   MODEL 

RH(%) 

PREDICTED 

ERR0R(%) 

_ 

OBS         ASSUMED 

3.  75|Um 

11.  00\im 

65              70 
75              70 
85              90 
95              90 

11.  8 

-10.  2 

8.  6 

-    5.  3 

9.  0 

-11.  8 

12.  6 

-13.  5 

b.     RURAL   MODEL 

RH(%) 

PREDICTED 

ERROR(%) 

OBS      ASSUMED 

3.  75fim 

11. OOum 

65              70 
75              70 
85              90 
95              90 

21.  7 
17.2 
13.  3 
28.  6 

16.  2 

23.  4 

9.  1 

30.  8 

E.       THREE-MODE  PARTICLE  SIZE  DISTRIBUTION  ERROR  ANALYSIS 

The  model  for  the  estimation  of  extinction  coefficients  assumes  that  the  marine 
particle  distribution  arises  from  particles  of  a  continental  and  marine  origin.  In 
conditions  of  high  winds  this  is  altered  by  the  introduction  of  larger  size  particles  into 
the  marine  boundary  layer.  As  the  wind  speed  increases  above  9  m/s,  waves  crest  and 
particles  are  ejected  into  the  atmosphere  (Monahan  er  al.,  1983).  These  particles  have 
an  average  radius  of  2  Jim  and  comprise  the  third  mode  of  the  particle  size  distribution 
(Gathman,  1983).  The  purpose  of  this  test  was  to  examine  the  effect  of  input  of  an 
extinction  coefficient  obtained  from  a  variable  wind  regime  and  estimate  extinction 
values  using  the  assumed  two-mode  model.  The  model  for  the  estimation  of  extinction 
coefficient  was  modified  to  allow  addition  of  a  third  wind-dependent  mode  represented 
by  variable  A3  as  shown  previously  in  Fig.  2.4.  The  values  of  A3  equal  to  1.0,  0.5  and 
0.1  were  taken  to  represent  the  production  of  10  Jim  particles  at  wind  speeds  of 
approximately  31.  27,  and  15  m/s  respectively.  The  shape  of  A3  is  described  by:  A3  = 
10exp(.06  x  wind  speed  -  2.8)  and  is  a  modification  of  the  Navy  Aerosol  Model 
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Fie.  4.3  Effect  of  changing  r2  on  the  distribution  of  marine  particles  by  ±20°o. 
Vertical  lines  at  Qcxt  peak  Tor  K  =  0.63  and  0.S6  ftm.  Sum  of  modes  1  &  2  (solid), 
mode  1  (dashed),  mode  2  (dot). 
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TABLE  6 

PARTICLE  SIZE  DISTRIBUTION  SENSITIVITY  STUDY 

MARINE 

MODEL 

RH(%) 

R2(%) 

PREDICTED 

ERROR(%) 

3. 

75nm 

11. 

OO^im 

50 
50 

50 

±    5 
±10 
±20 

0.  5 

-    4.  1 
-10.  5 

8.  7 
18.  2 
25.  0 

0.  0 
-    9.  7 
-20.  6 

11.  8 
24.  5 
40.  5 

70 
70 
70 

±    5 
±10 
±20 

0.  7 
-    4.  0 
-10.0 

8.  4 
17.  5 
22.  7 

-2.1 
-    9.  0 
-19.0 

11.0 
22.  9 
35.  8 

85 
85 
85 

±    5 
±10 
±20 

1.  3 
-1.1 
-8.2 

6.  4 
13.  6 
10.  7 

-1.3 
-    6.  1 
-17.  1 

9.  2 

19.  4 
14.  0 

90 
90 
90 

±    5 
±10 
±20 

-    5.  0 

8.  8 

14.  3 
21.  5 

-    4.  8 

11.  9 
21.  0 
36.  5 

(Gathman,  1983). [  The  model  was  run  with  TN  set  at  4,000  cm°  and  N2  set  at  0.02 
for  relative  humidity  values  of  70  and  90%  and  TN  was  set  at  4,135  cm  and  N2  at 
0.0135  for  relative  humidity  values  of  60  and  85%.  Different  parameters  were  selected 
to  test  the  sensitivity  when  matrices  were  not  interpolated.  The  procedure  was  as 
follows: 

1.  Obtain  an  initial  extinction  value  of  extinction  for  a  eiven  TN,  N2,  RH  and  A3 
value.    For  the  first  run  set  A3  to  1.0. 

2.  Enter  the  extinction  coefficients  into  the  first  section  of  the  prediction  model  to 
retrieve  a  value  for  TN  and  N2. 

3.  Enter  TN,  N2,  RH  and  an  A3  value  of  0.0,  which  reverts  the  program  back  to  a 
two-mode  distribution,  into  the  second  portion  of  the  prediction  model  to 
obtain  the  predicted  extinction  values  at  3.75  and  11.00  Jim. 

4.  Calculate  the  error  as  discussed  in  IV.  A. 

5.  Repeat  for  A3  values  of  0.5  and  0.1.  Also  repeat  for  several  relative  humiditv 
values  for  the  marine  model  only. 


From  H.  Hughes  by  personal  communication,  present  affiliation:  Ocean  and 

Atmospheric    Sciences    Division,    Naval    Ocean    Systems    Center,    San  Diego,    CA 

92152o000.                                                                    y                   .  B 
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With  the  addition  of  a  third  mode  we  observe  in  Table  "<  that  the  error  *s  greatest 
at  the  higher  wind  speed,  up  to  82%  for  A3=  1  and  \  —  11.00  Jim  and  up  to  39%  at 

a  =  3.75  fim.  Results  show  that  the  errors  decreased  as  relative  humidity  values 
increase.  The  addition  of  more  particles  at  the  higher  radii  for  lower  RH  values 
influences  the  result  more  and  thus  the  higher  error  at  the  lower  RH  values.  Also, 
since  the  higher  radii  are  weighted  more  by  the  extinction  efficiency  function  for  the 
higher  wavelength  a  higher  error  is  observed  at  11.00  \am.  There  is  a  slight  increase  in 
error  at  60%  RH  due  to  the  added  interpolation-  step  over  a  20%  range  of  relative 
humidity.  The  interpolation  error  is  not  observed  at  85%  RH  since  the  interpolation  is 
over  a  10%  spread  and  the  error  associated  with  the  introduction  of  the  third  mode 
dominates. 


TABLE  7 

EFFECT  OF  A  THREE-MODE  PARTICLE  SIZE  DISTRIBUTION 

MARINE 

MODEL 

RH(%) 

A3 

PREDICTED 

ERROR(%) 

3. 75nm 

11.  OO^im 

60 
60 
60 

1.  0 
0.  5 
0.  1 

-39.  0 
-30.  2 
-7.1 

-81.  7 
-75.  8 
-47.  2 

70 
70 
70 

1.  0 
0.  5 
0.  1 

-30.  2 
-22.  3 
-    4.  0 

-75.  5 

-67.  0 
-33.  1 

85 

85 
85 

1.  0 
0.  5 

0.  1 

-1.3 
8.  0 

-32*78 
-    3.  5 

90 
90 
90 

1.  0 
0.  5 
0.  1 

-  3.  7 

-  1.  3 

-  0.  1 

-34.  5 
-22.  4 
-    5.  9 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

A  method  for  estimating  extinction  coefficients  in  the  near-infrared  and  infrared 
wavelengths  using  channels  1  and  2  of  the  NOAA-7  AVHRR  was  developed.  The 
model  consists  of  two  sections.  Input  of  the  type  of  aerosol  model  (marine  or  rural) 
and  the  effective  relative  humidity  value  for  the  boundary  layer  defines  the  average 
mode  radius.  Input  of  the  extinction  coefficient  at  0.63  and  0.86  urn  determine  the 
magnitude  of  the  retrieved  parameters.  TN  and  N2,  necessary  to  describe  the 
distribution  of  aerosols.  The  second  section  uses  these  values  and  the  effective  relative 
humidity  value  for  the  boundary  layer  to  estimate  the  extinction  coefficients  at  3.75 
and  11.00  urn. 

Five  tests  were  developed  to  determine  the  limits  and  sensitivity  of  the  model. 
The  first  test  examined  the  error  associated  with  the  retrieval  of  a  'perfect'  satellite 
estimate  of  extinction.  Given  an  initial  value  of  extinction  at  0.63  and  0.86  Jim  for  a 
given  TN  and  N2,  results  show  that  the  retrieval  error  is  the  smallest  at  the  middle 
humidity  values  where  the  extinction  coefficient,  as  a  function  of  relative  humidity, 
changes  the  least  and  interpolation  occurs  over  a  relatively  small  range  (10%  spread). 

Results  for  the  second  test  show  that  the  models'  accuracy  is  basically  determined 
from  the  accuracy  of  the  input  extinction  coefficient  values.  Error  in  the  predicted 
extinction  is  less  than  3%  after  considering  the  error  in  input  extinction  for  all 
humidities  tested  and  for  both  models. 

Relative  humidity  affects  the  distribution  of  large  radii  particles  in  the  marine 
atmospheric  boundary  layer.  A  5%  error  in  the  value  of  the  relative  humidity  for  the 
area  of  interest  will  lead  to  errors  as  high  as  14%  for  the  marine  model  and  31%  for 
the  rural  model  in  the  estimated  extinction  coefficient  values  at  the  longer  wavelengths. 
This  test  shows  that  errors  in  relative  humidity  affect  the  rural  model  the  most,  since 
the  information  contained  in  the  second  mode  is  more  limited  for  the  rural  model  for  a 
given  relative  humidity. 

The  particle  size  distribution  can  also  be  changed  by  increasing  the  average  size 
of  the  particles  in  the  boundary  layer.  Thus  changes  in  the  mode  radius  by  5,  10,  20% 
will  result  in  an  error  in  the  estimated  extinction  coefficient  value  of  up  to  12,  24  and 
40%  respectively.  The  highest  errors  are  for  calculations  at  the  11.00  Jim  wavelengths 
since  the  extinction  efficiency  function  weighs  the -larger  radii  the  most. 
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The  last  test  evaluates  the  performance  of  a  two-mode  particle  size  distribution 
model  against  the  Nravy  aerosol  three-mode  model.  With  the  addition  of  a  third  mode, 
we  observe  that  the  error  decreases  as  the  wind  speed  decreases.  Since  the  third  mode 
contains  very  large  particles,  the  error  is  highest  for  estimates  at  11.00.  Jim  as  expected. 

Three  general  conclusions  can  be  made  from  this  study.  First,  estimates  of 
extinction  at  infrared  wavelengths  rely  heavily  on  the  accuracy  of  the  input  extinction 
coefficient  values  at  0.63  and  0.86  Jim.  Therefore,  care  must  be  exercised  in  obtaining 
the  best  possible  value  of  extinction  from  satellite  estimates  of  optical  depth  and 
NOLAPS  estimate  of  boundary  layer  depth.  The  second  is  that  extinction  is  dependent 
on  particle  size  distribution  and  thus  the  shape  of  the  distribution  used  in  the  model 
will  affect  the  estimated  value  directly.  And  finally,  factors  affecting  the  particle  size 
distribution,  such  as  relative  humidity  and  wind  speed,  must  be  taken  into  account. 

It  must  be  noted  here  that  the  model  as  it  stands  now  is  limited  to  certain  TN 
and  N2  values.  To  be  used  operationally  these  parameters  and  the  associated 
extinction  matrices  would  have  to  be  expanded.  Error  in  the  interpolating  technique 
can  be  decreased  by  increasing  the  number  of  relative  humidity  values  that  are  used  in 
the  model.  Also,  due  to  large  particles  being  generated  by  high  winds,  the  model 
should  be  modified  to  allow  for  variable  winds  similar  to  the  Navy  aerosol  model. 
Lastly,  methods  for  estimation  of  relative  humidity,  to  accuracies  of  95%  or  better, 
need  to  be  developed  for  meteorological  satellites  to  facilitate  measurement  of  relative 
humidity  in  remote  areas  and  thus  provide  better  input  into  the  model. 
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